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ABSTRACT 

The frequency of microlensing planet detections, particularly in difficult-to- 
model high-magnification events, is increasing. Their analysis can require tens of 
thousands of processor hours or more, primarily because of the high density and 
high precision of measurements whose modeling requires time-consuming finite- 
source calculations. I show that a large fraction of these measurements, those that 
lie at least one source diameter from a caustic or the extension from a cusp, can be 
modeled using a very simple hexadecapole approximation, which is one to several 
orders of magnitude faster than full-fiedged finite-source calculations. Moreover, 
by restricting the regions that actually require finite-source calculations to a few 
isolated "caustic features", the hexadecapole approximation will, for the first 
time, permit the powerful "magnification-map" approach to be applied to events 
for which the planet's orbital motion is important. 

Subject headings: gravitational lensing - planetary systems - methods numerical 



Introduction 



Microlensin g planets are bei ng discovered at an accelerating rate, with one being re- 



ported in 2003 ( 


Bond et al. 


2004 


). three in 2005 fUdalski et al. 


2005 


Gould et al. 


2006), two in 2006 


Gaudi et al. 


2008h. and perhaps as 



It has been a huge challenge for modelers to keep up with these discoveries, in large part 
because the computing requirements are often daunting: the parameter space is large, the 
surface is complex and generally contains multiple minima, and the magnification cal- 
culations are computationally intensive. Proper treatment of individual events can require 
tens of thousands processor hours, or more. Indeed, some potential planetary events have 
still not been fully modeled because of computational challenges. 



Planetary microlensing events are recognized through two broad channels, one in which 
the lightcurve perturbation is generated by the so-called "planetary caustic" that is directly 
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associated with the planet (iGould fc Loeblll992l ) and the o ther in which it is genera ted by 



the "central caustic" that is associated with the host star (iGriest fc Safizadehlll998l ). The 
former events are relatively easy to analyze and indeed the event parameters can be estimated 
reasonably well by inspection of the lightcurve. The latter events are generally much more 
difficult. 

There are several interrelated reasons for this. First, planets anywhere in the system 
can perturb the central caustic. This is why these events are avidly monitored, but by 
the same token it is often not obvious without an exhaustive search which planetary ge- 
ometry or geometries are responsible for the perturbation. Second, this very fact implies 
that several memb ers of a multi-planet system can be detected in central-caustic events 



(IGaudi et al.lll998l ). Multiple planets create a larger, more complicated parameter space, 
which can increase the computation time by a large factor. Even if there are no obvious 
perturbations caused by a second planet, an exhaustive search should be conducted to at 
least place upper limits on their presence. Third, if the source probes the central caustic, it 
is ipso facto highly magnified. Such events are brighter and more intensively monitored than 
typical events and so have more and higher- quality data. While such excellent data are of 
course a boon to planet searches, they also require more and more-accurate computations, 
which requires more computing time. Fourth, central-caustic planetary events are basically 
detectable in proportion to the size of the caustic, which roughly scales oc q/\b — 1|, where 
q is the planet/star mass ratio and h is the planet-star separation in units of the Einstein 
ring. Thus, these events are heavily biased toward planets with characteristics that make 
the caustic big. Such big caustics can undergo subtle changes as the planet orbits its host, 
and in principle these effects can be measured, thus constraining the planet's properties. 
Exploration of these subtle variations requires substantial additional computing time. Fi- 
nally, there is also a bias toward long events, simply because these unfold more slowly and so 
increase the chance that they will be recognized in time to monitor them intensively over the 
peak. Such long events often display lightcurve distortions in their wings due to the Earth's 
orbital motion which, if measured, can further constrain the planet properties. However, 
probing this effect (called "microlens parallax" ) requires yet another expansion of parameter 
space. Moreover, the "parallax" signal must be distinguished from "xallarap" (effects of the 
source orbiting a companion), whose description requires a yet larger expansion of parameter 
space. 

There are two broad classes of binary-lens (or triple-lens) magnification calculations: 
point-source and finite-source. The former can be used whenever the magnification is essen- 
tially constant (or more precisely, well-characterized by a linear gradient) over face of the 
source, while the latter must be used when this condition fails. Point-source magnifications 
can be derived from the solution of a 5th (or 10th) order complex polynomial equation and 
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are computationally very fast. When the point-source is well separated from the caustics 
(as it must be to satisfy the linear-gradient condition) then this calculation is also extremely 
robust and accurate. 

The main computational challenge in modeling planetary events comes from the finite- 
source calculations. Almost all integration schemes use inverse-ray shooting, which avoids 
all the pathologies of the caustics by performing an integration over the image plane (where 
the rays behave smoothly) and simply asks which of the rays land on the source. The 
problem is that a large number of rays must be "shot" to obtain an accurate estimate 
of the magnification, which implies that high-quality data demand proportionately longer 
computations. Of course, the higher the magnification, the larger the images, and so the 
more rays are required. There are various schemes to expedite inverse ray-shooting, including 
clever algorithms for identifying the regions that must be "shot" and pixelation of the source 
plane. The bottom line is, however, that the overwhelming majority of computation time is 
spent on finite-source calculations. 

Here I present a third class of binary-lens (or triple-lens) computation t hat is intermedi- 



ate be tween the two classes just described: the hexadecapole approximation. iPejcha fc Heyrovsky 



(120071 ) expand finite-source magnification to hexadecapole order and illustrate that the 
quadrupole term by itself can give quite satisfactory numerical results. In this paper, I 
develop a simple prescription for evaluating this expansion. While this algorithm is about 
10 times slower than point-source calculations, it is one to several orders of magnitude faster 
than finite-source calculations. The method can be applied whenever the source center is 
at least two source radii from a caustic or the extension from a cusp. Typically, well over 
half the non-point-source points satisfy this condition, meaning that the method can reduce 
computation times by a factor of several. Moreover, by isolating the small regions of the 
lightcurve where finite-source calculations must be used, the method opens up the possibility 
that the finite-source computations themselves can be radically expedited for the special, but 
very interesting, class of events in which planetary orbital motion is measured. 



2. Hexadecapole Approximation 

If the source does not straddle a caustic, then the magnification field can be Taylor 
expanded around the source center (xo,?/o) as a function of coordinate position {x,y) (all 
distances being expressed in units of the Einstein radius), 

oo n 
n=0 i=0 



-4- 



We wish to evaluate Afinito(p), the magnification of a source of radius p, 



J^dw J^"" dr]A{w,r]) 



(2) 



in terms of the Ani. Here {w,ri) are polar coordinates: (w cos?], w sin ?]) = {x — XQ,y — yo) 
To do so, we first average A{x, y) over a ring of radius w and obtain. 



where 



A{w) 



Aq = A, 



dr] A{w cos T], w sin rj) 
A20 + A22 



00, 



Ao + A2W'^ + A4W^ + . . . 

3^40 + Ai2 + 3A44 



A. = 



8 



(3) 



(4) 



To evaluate Afinite in terms of Aq, v42, ^4 . . ., we must first specify a limb-darkening law 
for the surface brightness S{w), which for simplicity we take to be linear. 



S{w) 



1-r l-:rWl-^ 



(5) 



where F is the limb- darkening coefficient and F is the source flux. We then substitute into 
equation ([3]) to obtain 



finite 



j,Uw2nwA{w)S{w) _ A^pV 1 X A,p' ^ 11 

F ° 2V 5/ 3V 35' 



Note that the "limb-darken ing factors" in parentheses are simply {2ti / F) dwS{ 



w jw 



(6) 



2n+l 



(jPejcha fc Heyrovskyll2007l ). and therefore do not depend on the magnification in any way. 
Hence, for any given adopted limb-darkening profile, these can be calculated just once. 

Let us now assume that the field is adequately described by a hexadecapole. Averaging 
over the four points on a ly-ring that are shifted by an arbitrary angle relative the cardinal 
directions, we obtain 



1 ^ 

j=0 



w sm 



TT 

^2 



A, 



A2W^ + 



(v44o + v444) (1 + cos^ 20) + v442 slu^ 20 



(7) 



while rotating this geometry by 7r/4 gives. 



2 (A4O + A44)(l + sin2 20) + A42COs2 20 ^ 

^«),x = A2W H ^ w . 
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For a given source of size p and position {xq^uq), one can therefore determine Aq, Ap,+, 
5 and Ap,x5 from a total of 13 point-source calculations, and thus derive 

= ~ = i-±±^ - A,p\ (9) 

which can then be substituted into equation IQ to obtain Aa^iteip, xo,yo). 



3. Range of Validity 

Clearly, this approximation cannot be used when the source lies on a caustic, but how 
close can the source be before the approximation breaks down? The breakdown will be 
driven by the leading term of the caustic's singularity, so it is sufficient to examine idealized 
cases whose leading-order behavior is the same as that of real caustics. More formally, one 
can write the magnification field as the linear sum of an idealized singularity and a more 
complicated, but well-behaved field. Only the former will contribute to the breakdown of the 
approximation. Since equations (JH]), dZ]), dHD, and are strictly linear, this decomposition 
is absolutely rigorous. Caustic singularities come in two basic varieties, fold caustics and 
cusps. The former diverge as (Am_l)~^/^ as one approaches the caustic from the inside, where 
Au± is the perpendicular distance to the fold. The latter are much more complex. They 
diverge as {Au)~^ in the immediate neighborhood of the cusp as one appr oaches it from 



the outside. However, at greater distances, they develop into long "fingers" (ICould fc Loeb 
I992I : IPeicha fc HevrovskiboOTh . 



3.1. Fold Caustic 

For simplicity, I consider a uniform source that lies a distance zp (where p is the source 
radius) from the fold. I begin by assuming that the magnification is dominated by the two 
"new images" that meet on the critical curve. I then discuss how the result is changed when 
this assumption is relaxed. The magnified fiux is then given by 



V "J- m=\} 

which after some algebra [e.g., dyy'^il — y)^ = m\n\/{m + n + 1)!], simplifies to 



_- (4n- l)!!(2n -!)!! _,, 

^ 23"(2n)!(n + l)! ' ^ ' 



where n = m/2. That is. 



Aiz) 3 2 35 4 1155 45045 o , , 
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Hence the error due to the hexadecapole approximation is 

AA{z) 1155 _g 



2.8x10-^(0" . (13) 



If one is "sufficiently near" a fold caustic, the magnification will always be dominated 
by the two "new images". However, for planetary caustics, this will generally no longer be 
the case even at one or two source radii from the caustic. The net effect is to "dilute" the 
caustic divergence and so to make the the hexadecapole approximation valid at even closer 
separations than indicated by equation ( fT3l) . See §[5l 

3.2. Cusp 

For simplicity, I analyze the case of a {Au)~^ divergence and then discuss how the results 
may be expected to change for real cusps. In this case, we have 

m=0 1=0 

which after some algebra [e.g., ((2 cos 6')^''') = C|'^] reduces to 

^ ^ - ,~2n n 1 (2n + 2i-l)!! 

where, again, n = m/2. That is, 

A{z) ^ 1 2 3 4 25 6 245 o 

so that the error due to the hexadecapole approximation is, 

AA(z) 25 . Afz\-^ 

— = 3.8 X 10-^ - . (17) 



Ao 210 V2 

Equation (fT71) gives a reasonable lower limit on how closely one may approach a cusp 
using the hexadecapole approximation. However, because the cusp develops a linear, finger- 
like structure at moderate distances, this approximation can fail well away from the cusp. 
Nevertheless, I find numerically that the total duration of the failure is generally represented 
reasonably well by equation (fT7|) . See §0 

Thus, for both fold caustics and cusps, the hexadecapole approximation can reduce the 
error below 0.1%, except for intervals characterized by z = 2.5. 
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4. Implementation 

The main consideration wlien implementing this approximation is to make certain that it 
is apphed only in its range of validity. This is easiest when one is probing an already-located 
minimum, which is often the most time-consuming part of the investigation. One can then 
simply plot the difference between the hexadecapole approximation and the finite-source 
calculation for a single model as a function of time, and so locate empirically the regions of 
the former's range of validity. If a "safety zone" is placed around the finite-source regions, 
then there is little danger that it will be crossed during the minor excursions that occur 
while probing a minimum. It is straightforward to determine automatically whether such 
unexpected crossings are occurring simply by comparing the finite-source and hexadecapole 
calculations for the first and last points of each finite-source region. If these do not agree, 
the "safety zone" has been crossed. 

One must be more careful when applying this method to blind searches because it is 
harder to determine whether any particular stretch of the lightcurve is either crossing or 
nearby a caustic. In some cases, this will be straightforward and in others more difficult. 
The one general point to note is that the same "safety zone" check can be made. 

Finally, I remark that in some of the lightcurve regions where the hexadecapole ap- 
proximation is valid, it may be overkill: the quadrupole or even monopole (point-source) 
approximations may be perfectly satisfactory. Again, it is straightforward to find these 
subregions, simply by mapping the hexadecapole/quadrupole and hexadecapole/monopole 
differences for a single model. 

5. An Illustration of the Method 



Figures [T] and [2] give a practical example of the hexadecapole approximation. The top 
panel of Figure [T] shows a caustic geometr y (black) that is ba. sed loosely on the caustic in 



the planetary event OGLE-2005-BLG-071 (jUdalski et al.l 120051 ). but with a different source 



trajectory (blue) that has been chosen the maximize the number of illustrative "features": 
the source passes by two cusps and then enters and exits the caustic. 

The middle panel shows the resulting lightcurve (black) together with three successive 
levels of approximation: monopole (i.e., point-source blue), quadrupole (red), and hexade- 
capole (green). The bottom panel shows the residuals of the latter three relative to the 
first. 

There are several points to note. First, the hexadecapole approximation works ex- 
tremely well over the entire region shown except for a few source-radius crossing times in 
the immediate vicinity of the first cusp approach and the two caustic crossings. Second, the 
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point-source approximation basically does not work at all over the entire region shown, at 
least if one is attempting to achieve precisions of 0.1% (which is typically required). Third, 
there are significant regions where the quadrupole approximation is adequate. 

Figure [2] is a zoom of Figure [1], focusing on the region of the two caustic crossings. The 
crosses in the bottom panel show the predictions for the breakdown of the quadrupole (red) 
and hexadecapole (green) approximation inside the caustic (eqs. [12] and [I3]), assuming 
one is trying to achieve a precision of 0.05%. The black crosses indicate the prediction for 
breakdown outside the caustic, namely one source radius from the caustic. 

Both the quadrupole and hexadecapole approximations prove to be too conservative in 
that the approximations remain valid substantially closer to the caustic than expected. The 
reason for this is clear from the middle panel: at the points of the predicted breakdown, the 
underlying magnification profile is no longer dominated by the square-root singularity that 
is produced by the two "new images" . Recall from § 13.11 that equations (fT2l) and (fT3|) were 
specifically derived under the assumption these images would dominate the magnification. 

Returning to Figure [H the situation is more complex for the cusp approaches. Equation 
f fT7|) predicts that the hexadecapole approximation will break down from —25.8 to —22.7. 
The actual breakdown is from —27.6 to —23.2. This displacement to the left reflects the 
"finger" of enhanced magnification that extends from the cusp axis, as indicated in the top 
panel. In fact, I find that for trajectories passing farther from this caustic, the breakdown 
occurs at even earlier times, again corresponding to the intersection of the trajectory with 
the cusp axis (rather than the point of closest approach). Moreover, the breakdown continues 
to occur even when the point of closest approach is well beyond the separation predicted 
by equation (ITTI) . Nevertheless, the full duration of the breakdown is never much greater 
than 2z as predicted by this equation. In brief, for cusps, the hexadecapole approximation 
does work except for brief intervals whose duration is given by equation (IT7|) . but the time 
of breakdown is not accurately predicted by this equation. Similar remarks apply to the 
quadrupole approximation. 

6. Application to Lenses With Orbital Motion 

While the hexadecapole approximation can save substantial computation time in a wide 
range of cases, it may be especially useful for planetary lenses with meas urable orbital mo - 



tion by rendering them accessible to the "magnification map" technique (iDong et al.ll2006l ). 
Magnification maps are potentially very powerful. One first chooses "map parameters", i.e., 
(6, q) for a single planet or (fei, gi, 62, Q'2, 4>) for a triple system, where is the angle between 
the two planets. Then one shoots the entire Einstein ring (out to a specified width corre- 
sponding to, say, magnification A = 100). One then both stores the individual rays and 
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tiles the source plane with hexagonal pixels, keeping track of the number of rays landing in 
each pixel. Pixels landing wholly within the source are evaluated at their centroid, while 
pixels that cross the boundary are evaluated on a ray-by-ray basis. Using this map, one can 
minimize over the remaining microlensing variables (time of closest approach to, impact pa- 
rameter Mo, Einstein timescale ^e, source trajectory angle a, source size p, as well as parallax 
and xallarap parameters, if these are needed). Each map can be created in a few seconds 
and fully explored in a few minutes, thereby permitting a rapid Markov-chain walk through 
map space. 

The drawback is that, to date, magnification maps have not been applicable to lenses 
with significant orbital motion: to the extent that the lens separation changes during the 
event, different maps would be needed at different phases of the event, potentially a very 
large number of them. In some cases, this problem is now completely resolv ed. For example. 



OGLE-2005-BLG-071 exhibits some signature of rotation (IDong et al.ll2008l ). which had been 
difficult to probe simultaneously with finite-source effects. Because the source trajectory 
comes no closer than z = 10 source radii from the cusp, this event can now be handled 
completely in the hexadecapole approximation. 

However, even for events with one or several caustic features that require finite-source 
calculations, it will now be possible to evaluate these with maps. Each caustic feature lasts 
only about the time required for the source to cross its own diameter, which is typically 
a few hours. The orbital motion during these features is negligible, implying that the lens 
geometry can be adequately represented by a single map. Several maps can be created, one 
for each feature occurring at different times. The points "between features" can be evaluated 
in the hexadecapole approximation, which allows a continuously evolving lens geometry. 

7. Conclusions 



I have identified an intermediate regime between the ones where finite-source effects are 
dominant and negligible, respectively. In this regime, magnifications can be evaluated with 
very high precision using a simple hexadecapole approximation, for which I give a specific 
prescription. Outside of small (few source-diameter crossing times) regions associated with 
caustic crossings and cusp approaches, the approximation has a fractional error of well under 
0.1%. Some events can now be analyzed without any traditional finite-source calculations, 
while for others these calculations will be drastically reduced. In particular, by restricting 
the regions that do absolutely require finite-source calculation to a few isolated zones, this 
approximation opens the possibility of applying "magnification maps" to planetary systems 
experiencing significant orbital motion. 
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Fig. 1. — Top: source trajectory (blue) through caustic geometry (black) of a simulated high- 
magnification microlensing event, in units of the source radius. Middle: resulting lightcurve 
as found from full finite-source calculation (black) and the monopole (blue), quadrupole (red), 
and hexadecapole [green) approximations. Bottom: residuals of the three approximations 
relative to the full calculation. See § |5] for full discussion. 
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Fig. 2. — Zoom of Fig. [T], focusing on region of the two caustic crossings. Crosses in middle 
panel show predictions for range of validity of the quadrupole {red) and hexadecapole (green) 
approximations inside the caustic and for both (black) outside. See §[5] for full discussion. 



